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Abstract: We study the evolution of entanglement entropy in a 2-dimensional equilibra- 
tion process that has a holographic description in terms of a Vaidya geometry. It models 
a unitary evolution in which the field theory starts in a pure state, its vacuum, and un- 
dergoes a perturbation that brings it far from equilibrium. The entanglement entropy in 
this set up provides a measurement of the quantum entanglement in the system. Using 
holographic techniques we recover the same result obtained before from the study of pro- 
cesses triggered by a sudden change in a parameter of the hamiltonian, known as quantum 
quenches. Namely, entanglement in 2-dimensional conformal field theories propagates with 
velocity v 2 = 1 [1]. Both in quantum quenches and in the Vaidya model equilibration is 
only achieved at the local level. Remarkably, the holographic derivation of this last fact 
requires information from behind the apparent horizon generated in the process of gravita- 
tional collapse described by the Vaidya geometry. In the early stages of the evolution the 
apparent horizon seems however to play no relevant role with regard to the entanglement 
entropy. We speculate on the possibility of deriving a thermalization time for occupation 
numbers from our analysis. 
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1. Introduction 

The AdS/CFT correspondence has proven to be a unique tool for the study of strongly 
coupled quantum field theories. Since its first proposal a fruitful subject of application 
has been thermodynamics, where the holographic dictionary relates the plasma phase of 
strongly coupled field theories at thermal equilibrium to black hole geometries [2]. In 
the last years the correspondence has been generalized to describe small deviations out of 
equilibrium in the context of linear response theory, and especially its hydrodynamic long 
wavelength limit [3]. It turns out that holographic plasmas are close to perfect fluids, with 
very low specific viscosity [4]. This fact describes well the data from heavy ion experiments 
and has represented a major success for the holographic techniques, which were extended 
in [5,6] to the study of non-linear fluid dynamics. 

Continuing this line of work, it is very interesting to apply holography to the analysis of 
strongly coupled field theories far from equilibrium. On general grounds, a system brought 
into a highly excited state is expected to evolve towards a stationary state that at the 
macroscopical level can be described as thermal equilibrium. The holographic counterpart 
of this process should correspond to a process of gravitational collapse ending in the for- 
mation of a black hole [7-9] . Recently, the far from equilibrium evolution of an anisotropic 
plasma [10] and a boost invariant plasma [11, 12] have been addressed using holography. 
In [13] the collapse of a shell build out of an AdS massless scalar was analyzed in the limit 
of weak fields, discussing its implications for field theory thermalization at strong coupling. 
Time-dependent backgrounds associated with a localized sector of the dual field theory 
have been constructed in the probe approximation in [14]. Articles [10-13] mostly consid- 
ered field theory observables that can be extracted from an expansion of the dual gravity 
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background close to the boundary, in particular the energy-momentum tensor. We want 
instead to follow the thermalization process using a quantity for which this information is 
not enough. 

Entanglement is a fundamental notion in quantum physics. A measure of quantum 
entangled in extended systems is provided by entanglement entropy. The reduced density 
matrix of a subsystem A is obtained by tracing over the degrees of freedom of its comple- 
mentary, pa = Tr^ p, where p is the density matrix of the entire system. The entanglement 
entropy of region A is defined as the von Neumann entropy of this reduced density matrix 



When the system is in a pure state, p = it gives a measure of the quantum 

entanglement between the degrees of freedom inside and outside A. In this situation the 
von Neumann entropy of the complete system is zero and the following property is fulfilled 



When the system is instead in a mixed state, the entanglement entropy receives also sta- 
tistical contributions. 

The study of the entanglement entropy is a subject of active research. Many results 
are available for 2-dimensional field theories and quantum spin chains, see [15] for a re- 
view and references. Dynamical situations have been also considered, in particular that 
triggered by a sudden change in a parameter of the hamiltonian [1]. This action, known as 
quantum quench, can be realized experimentally in systems of cold atoms, together with 
a subsequent time development under negligible dissipation. This has provided a strong 
motivation for the study of quantum quenches. In this set up, the system is prepared in 
a pure state before the quench, |\I>o)> and due to the unitary character of the subsequent 
evolution it remains in a pure state, | As mentioned above, under this condition the 
entanglement entropy will describe the evolution of quantum entanglement in the system. 
General results were obtained in [1] for 2-dimensional conformal field theories. It was shown 
that entanglement propagates from small to large scales with velocity v 2 = l. This has the 
important consequence that equilibration can only be achieved at the local level. 

We will address the evolution of entanglement entropy using holographic techniques. 
The AdS/CFT correspondence has proven to be an efficient method for the calculation of 
the entanglement entropy in those field theories that admit a gravity dual, see [16] for a 
review. The initial proposal for the holographic derivation of entanglement entropy in static 
backgrounds [17,18], based on the calculation of minimal surfaces, was generalized to time 
dependent situations in [19]. We will focus on the duality between 2-dimensional CFT's 
and gravity on asymptotically AdS3 geometries. In this case the holographic proposal 
relates the entanglement entropy with the length of spacelike geodesies that anchor on the 
AdS boundary. We will use a Vaidya geometry, which describes the formation of a black 
hole out of a collapsing shell of null dust, as dual model for a field theory equilibration 
process. The Vaidya metrics are analytic solutions of Einstein equations that will allow us 
to carry out a detailed analysis of the relevant geodesies. 



Sa = -Tr A (pa In pa) • 



(1.1) 



Sa = Sa ■ 



(1.2) 
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This holographic model realizes a similar set up as the quantum quench. The system 
starts in a pure state, in this case the vacuum of the dual CFT, and undergoes a pertur- 
bation, modeled by the collapsing shell, which brings it into a highly excited state. We 
will argue that (1.2) is verified along the complete evolution, implying that the evolution is 



unitary and the system always remains in a pure state. The Vaidya model and the quantum 
quench only differ in the type of perturbation that triggers the time development, which 
determines the initial entanglement pattern. While for the quenches considered in [1] en- 
tanglement is initially localized on very small scales, in the Vaidya model there are long 
range correlations in the early stages of the evolution. However we will recover the same 
conclusions about the evolution of entanglement in 2-dimensional CFT's. This offers a new 
and strong test of the AdS/CFT correspondence, and at the same time provides a new ex- 
ample where a detailed analysis of the time-development can be performed. A holographic 
model for a quantum quench that only affects some localized degrees of freedom has been 
proposed in [14]. 

An interesting question for the AdS/CFT correspondence is the possible role of gener- 
alized notions of horizon [20,21], like the apparent horizon, in the holographic description 
of field theory dynamics far from equilibrium. The Vaidya model will allow us to explore 
the relevance of the apparent horizon for the description of the dual field theory entangle- 
ment properties. Contrary to the case of the event horizon in static geometries, we will see 
that the apparent horizon does not act as a wall for the geodesies relevant to the derivation 
of the entanglement entropy. As happened for quantum quenches, our results show that 
relaxation is only achieved at the local level. Interestingly geodesies able to detect this 
effect reach behind both the event and the apparent horizons. Moreover, in the first stages 
of the evolution the geodesies do not appear to feel the location of the apparent horizon. 
Hence the apparent horizon seems to play a weak role for the description of quantum 
entanglement in the far from equilibrium field theory evolution. 

The organization of the paper is as follows. Section 2 contains a brief summary about 
entanglement entropy in 2-dimensional CFT's and its evolution after a quantum quench, 
together with the holographic proposal for its evaluation. In section 3 we introduce the 
Vaidya geometry. In section 4 we present the numerical study of spacelike geodesies that 
anchor on the boundary of a 3-dimensional AdS- Vaidya geometry, with emphasis on the 
role played by the apparent horizon. The implications for the evolution of entanglement 
are analyzed in section 5, comparing with the known results for quantum quenches. We 
end in section 6 with a further discussion, where we comment on the possibility of deriving 
a thermalization time for occupation numbers from our results. 

2. Entanglement entropy 

A remarkable property of entanglement entropy is its universal behavior in 2-dimensional 
conformal field theories [22-26]. In the ground state of a CFT living on a line, the entan- 
glement entropy of an interval scales logarithmically with its size 

S(l) = | log l -, (2.1) 
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where e is a UV cutoff. The only information relevant to S(l) is the CFT central charge, c, 
a quantity that provides a measurement of the number of degrees of freedom in the theory. 
In ( |2.1| ) we have omitted an /-independent non-universal contribution. Up to the same 
constant, the entanglement entropy of an interval at thermal equilibrium is given by the 
following expression [26] 

S,(0 = flog (iL.bhI*), (2.2) 

where f3 = 1/T is the inverse of the temperature. The ground state result is recovered 
for intervals much smaller than /3. For larger intervals the cutoff independent part of the 
entanglement entropy behaves as an extensive quantity and leads to the same result as the 
thermal entropy. 

The evolution of entanglement entropy in 2-dimensional systems after a quantum 
quench was analyzed in [1]. In this set up the system is prepared in its ground state, 
and at t = a parameter of the hamiltonian is suddenly changed. The old ground state 
provides a highly excited starting point for the evolution with respect to the new hamil- 
tonian. Explicit results can be obtained for 2-dimensional field theories when the ground 
state previous to the quench has a mass gap while the dynamics after the quench is con- 
formal. For large times and intervals, i, I » To with To of the order of the inverse mass gap 
in the initial state, the entanglement entropy of an interval is given by [1, 15,27] 



s(i,t) * Zh£± + { 



net , ,„ 
t<l/2. 



6r 

TTCI 



(2.3) 

t>l/2. 



12r 

The sharp cusp in this expression at t = l/2 should be understood as smoothed out over a 
region of order tq. After the quench S(l,t) grows linearly with time until it saturates at 
t~l/2. Comparison with (2.2) shows that for late times the entanglement entropy behaves 
as that of a thermal state, with 4to playing the role of the inverse temperature. An 
analogous conclusion is derived from the evaluation of correlation functions [27, 28] . This 
leads to interpret 1/4to as an effective temperature for the long wavelength modes [27]. 

The physical interpretation of ( |2.3[) is as follows. The non-zero energy density created 
at t = by the quench translates into the production of a translationally invariant sea of 
excitations. Entangled excitations will tend to spread over ever larger regions mainly as 
a result of propagation. Since the ground state before the quench has a mass gap and 
hence a finite correlation length, only excitations produced approximately at the same 
point will be entangled. A simple picture where entangled excitations are assumed to move 
classically and without scattering, was further proposed in [1]. For the critical quench 
described above, their velocity is taken to be v 2 = 1 equal for all momentum modes. At 
time t left and right-moving entangled excitations will have separated a distance 2t. In 
order to contribute to the entanglement entropy, only the left-moving or the right-moving 
components must lie inside the chosen interval, see Figjl|. When the separation reached 
by entangled excitations is smaller than the size of the interval, the number of them that 
contributes to S(l,t) grows with time proportional to At. When instead 2t>l no left-right 
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moving entangled excitations are contained inside the interval and S(l,t) saturates at a 
value proportional to 21. Hence, although causality prevents the system to equilibrate 
at the global level, any finite region will locally relax to a stationary state for which its 
complementary region acts as a thermal bath. 




Figure 1: Derivation of the entanglement entropy of an interval according to the model proposed 
in [1]. Only excitations originating from the shaded region at t — contribute to the entanglement 
entropy of region A at time t. 

The evolution of the entanglement entropy after a quantum quench has also been 
studied in spin chain models [1,29-32]. As in the previous case, it was found in several 
examples that the entanglement entropy of a block of spins grows linearly with time until 
it saturates at a value set by the number of spins in the block. However the way of 
approaching the saturation value shows that, while there is a maximal velocity for signal 
propagation, in general not all modes responsible for the propagation of entanglement reach 
it. The evolution ( |2.3D , modeled in Fig.^J, implies instead that in a 2-dimensional CFT the 
entanglement pattern of the excited state created by the quench is propagated by all modes 
with maximal velocity, v 2 = l. 

The evaluation of the entanglement entropy is a very complicated problem in general 
situations, since it involves performing traces of the density matrix. Therefore having 
alternative tools to address this problem is very valuable. This has been provided in the 
last years by the AdS/CFT correspondence for those field theories that admit a gravity 
dual [17,18]. The proposal referred to static situation, where there exists a timelike Killing 
vector in terms of which to define a foliation of the dual geometry by slices of constant 
time. Restricting to such a slice, the entanglement entropy of a region A is calculated 
holographically as the area of the minimal surface 7 that anchors at dA, the boundary of 
A, at the boundary of the dual gravity background 

Area(-f) 

Sa = ' {2A) 
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with Gm the bulk Newton constant. This proposal has been generalized to time dependent 
situations, where there is no unique way to extend a constant time slice of the boundary 
into the bulk geometry. It has been argued however that (|2.4| ) holds, interpreting 7 as 
an extremal surface in the unrestricted gravity background subject to the same boundary 
conditions [19]. 

In the absence of other scales than the temperature, the dual gravity description of 
a field theory in thermal equilibrium is given by a black hole geometry [2]. We want to 
study the dynamics of field theories far from equilibrium using holographic techniques. 
On general grounds the field theory will tend, as time evolves, to a state that at the 
macroscopical level can be described as thermal equilibrium. The candidate gravity duals 
should therefore represent a process of gravitational collapse ending in the formation of 
a black hole [7-9]. There is a family of exact solutions to the Einstein equations which 
describes the collapse of null dust to form a black hole, known as Vaidya metrics [33]. 
They were derived for asymptotically flat spaces, but can be generalized to spaces with a 
cosmological constant. In this paper we will focus on an equilibration process that can be 
described holographically in terms of a Vaidya geometry. 

More general collapse processes, triggered by an inf ailing shell built up from the non- 
normalizable mode of an AdS massless scalar were considered in [13]. According to the 
holographic dictionary, they describe the perturbation of the dual conformal field theory 
out of its vacuum state by switching on a time-dependent coupling to a marginal operator. 
The solution of Einstein equations in that case could not be obtained in a close form, being 
approached instead in terms of a perturbative expansion in the limit of weak perturbations. 
Since null dust models incoherent massless radiation, the Vaidya geometry can be expected 
to represent a limiting case of the previous set up. Constructing a holographic dual for a 
quantum quench would require a different dual gravity background, which includes an IR 
cutoff able to represent a mass gap in the field theory previous to the quench. Therefore the 
Vaidya geometry offers us the possibility to study the evolution of entanglement entropy 
in a different dynamical process to that generated by a quantum quench. 



3. Vaidya geometry 

The Vaidya metrics in an asymptotically AdS spacetime of dimension d+1, in Poincare 
coordinates are given by 

ds 2 = -(r 2 - ^\ dv 2 + 2drdv + r 2 £ dx 2 . (3.1) 
^ r ' i=l 

They are solutions of Einstein's equations in the presence of the following energy-momentum 
tensor 

m , si dm 

with all other components being equal to zero. The energy-momentum tensor characteristic 
of null dust satisfies oc k^k u with k 2 = 0. The previous expression is of this form with 
k/j, = 6/j, v , which generates infalling radial null geodesies. The metric (3.1) describes the 
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formation of a black hole out of collapsing null dust, or more in general, the absorption of 
null dust by a black hole. The Vaidya metrics were first derived as a model for the exterior 
geometry of a radiating star [34,35], being formulated in that case in terms of an outgoing 
null coordinate. 

In the study of gravitational collapse, a central question is how to characterize the 
formation of a black hole at the level of the local time evolution. The event horizon, 
being a global property of the spacetime, is not appropriate for this purpose. Several 
generalized notions of horizon have been proposed based on identifying black holes with 
regions containing trapped surfaces [20,21]. A trapped surface is a closed spacelike surface 
of co-dimension two for which the expansions along the two future directed normal null 
directions are negative. These generalized notions of horizon have in general the drawback 
of requiring a choice of spacelike surfaces and therefore not being covariantly defined. For 
the Vaidya metrics a natural choice are surfaces that respect the translational invariance 
along the X{ coordinates. With this choice most proposals coincide. One of them is the 
apparent horizon, defined as the boundary of trapped surfaces associated to a given foliation 
[36]. Under some smoothness conditions, the apparent horizon on each leaf of the foliation 
coincides with the co-dimension two surface for which the expansion of outgoing future 
directed normal null geodesies is zero, while that of ingoing ones contracts. A foliation 
preserving invariance under translations on the x% coordinates selects the surfaces r = const, 
v = const. 1 . Their associated normal null vectors are 

N in = -d r , N out = d v + -[r 2 - -^\d r (3.3) 

whose normalization has been chosen such that Ni n ■ N ou t = —1. The expansion along 
both null directions is given by 6 = g^V l _ l N u , where g^ v = g^ u + N^ n N^ ut + N^ ut N^ n is the 
induced metric. We obtain 

*. =— . *- = ^(r m 

and hence the apparent horizon is located at 

r h = {m{v)) 1 ' d . (3.5) 

We will also call apparent horizon the co-dimension one surface defined by the union of the 
surfaces 

Satisfying the null-energy condition [36], T^N^N" > for any null vector N^, implies 

dm 

->0. (3.6) 

Under this condition, the apparent horizon always grows. We will assume that the mass 
function m(v) tends to a finite value, m, as v—>oo. The location of the event horizon is 
given by 



1 These surfaces could be rendered closed by compactifying the x% coordinates. 
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It always lies at larger values of the radial coordinate than the apparent horizon, see Fig.|2|. 
Contrary to the apparent horizon, the event horizon is different from zero well before the 
collapse process starts. 




Figure 2: Position of the event (solid) and the apparent (dashed) horizons of the Vaidya metric 
~Jh for d = 2 and m(u) = (taaht>+l)/2. 



When the mass function is constant, m(v) = m, ( |3.l[) reproduces the Schwarzschild 
back hole geometry with planar horizon 

ds> = - (r 2 - dt 2 + + r 2 g d X \ , (3.8) 



i=i 



as can be seen by the coordinate change 



df r d ~ 2 

v = t + f(r), j- = - d . (3.9) 

dr r a — m 



The relation between the coordinates v and t in the (|3.8|) geometry, close to the boundary 
and the event horizon is 

I 

, r — > oo 

r 

v ~ t-\ 3_ i n ( r _ m yd\ ( r m i/d (3 1Q ) 

dm 1 '" 1 

Close to the boundary the foliations defined by both coordinates asymptotically coincide 
and define the same boundary time. However they differ strongly in the interior. In 
particular close to the event horizon t finite implies v— >— oo and conversely v finite t— >oo. 

In the following we will focus on the AdSs/CFT2 duality. For static backgrounds, like 
AdS3 and the BTZ black hole with non-compact horizon described by (3.8) with d=2, the 



holographic evaluation of the entanglement entropy requires to study spacelike geodesies 
that live in a constant t slice [17, 18]. The central segment of the geodesic will probe 
smaller values of the radial coordinate as the distance between its endpoints, anchored at 
the AdS boundary, is increased. Since v is a monotonously increasing function of r outside 
the event horizon, the geodesies will at the same time access earlier and earlier values of 
v. Moreover, those long enough to approach the event horizon, using ( 3.10| ), must reach 



-8- 



v— >— oo. We will show in the next section that geodesies evaluated in the Vaidya geometry 
have a similar behavior in terms of the v coordinate. Namely, for any fixed boundary time 
they explore arbitrarily early values of v. The map that this implies between measurements 
at a given time in the dual field theory and bulk physics will determine the evolution of 
the entanglement entropy. 



4. Spacelike geodesies in AdS 3 -Vaidya 

The AdS3-Vaidya metric with boundary R 1,1 is given by 

ds 2 = -(r 2 - m(v))dv 2 + 2drdv + r 2 dx 2 . (4.1) 

Since we are interested in describing a black hole formation process, m(v) must interpolate 
between and a finite value m, which sets the final black hole temperature T = y / m/27r. 
In this paper we will consider the analytic function 

, , tanhfv/a) + 1 ... 
m(v) = m ^-LJ. . (4.2) 

We will denote the boundary value of v by t, the field theory time. This geometry provides 
the holographic dual of a strongly coupled 2d CFT on a line in its ground state for t<—2a. 
Around t = the ground state is perturbed in a translationally invariant way creating a 
finite energy density that smoothly builds up in a time extent 5t ~ 4a. The system will 
then evolve towards a state of thermal equilibrium at the temperature set by the black hole 
collapse process. The expectation value of the field theory energy-momentum tensor can 
be read from the variation of the dual gravity action, properly regularized, with respect to 
the boundary metric [37, 38] . We obtain for the field theory energy density and pressure 



e(t) = P (t) = m(t) . (4.3) 

In [19] it was proposed that the holographic derivation of the entanglement entropy 
[17,18] extends to dynamical backgrounds. The entanglement entropy for an interval of size 
I is then given by the length of the spacelike geodesic, most conveniently parameterized by 
the coordinate x G [0, 1], and described by functions r(x) and v(x) subject to the boundary 
conditions 

r(0) = r(0 = oo, v(Q)=v(l)=t. (4.4) 
With this parameterization, the length functional for AdS3-Vaidya is 

L= dx\J r 2 + 2rV — (r 2 — m(v ))v' 2 . (4-5) 
Jo 

The absence of explicit x-dependence in the integrand leads to the following integral of 
motion 

r 4 

_ =r 2 + 2rV- (r 2 - m(v))v' 2 , (4.6) 
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where r* is a constant. The two second order differential equations obtained from extrem- 
izing ( f4,5| ) can be reduced with the help of the previous conservation equation to 



2 2/2 // 

r — r v — rv 



+ 2wV = . (4.7) 



The boundary conditions (4^4), together with the translational invariance of the metric in 
the x direction, imply that the associated geodesies will be symmetric with respect to the 
midpoint, at which r' = v' = Q. For this type of geodesies the integration constant r* gives 
the value of r attained at the midpoint. The previous differential equations with the mass 
function ( |4.2| ) can only be integrated numerically. Some results for the case of compact x 
were already presented in [19]. Without loss of generality, we will perform the numerical 
integration of the geodesic equations for a final black hole mass m = 1, or equivalently 
f5 = 2ir. The parameter m can be reintroduced by a simple rescaling, under which the time 
extent of the perturbation transforms as a^-a/y/m. 

Substituting ( |4.6| ) into the length functional and using the symmetry of the geodesies, 
we obtain the simplified expression 

2 r 1 / 2 

L(l,t) = — / r(x) 2 dx. (4.8) 



r* Jo 

The equations ( f4.6| ) , ( [4.7D for large values of the radial coordinate reduce to those describing 
geodesies in AdS3. As in that space, the leading behavior close to the boundary is given 
by 

r * 16 

X -^2< V ~ t ~- + ^2' ( 4 - 9 ) 

where we have focused on a neighborhood of the geodesic endpoint x = 0. The integration 
constants r* and 6 are to be fixed by requiring that the boundary conditions (4.4) are 



fulfilled at x = l. This implies 6 = both for AdS3 and the BTZ black hole, while for the 
dynamical Vaidya geometry we will show below that this is in general not the case. The 
behavior of r(x) close to the boundary gives rise to a logarithmic divergence of the geodesic 
length, which can be regularized by cutting the integration in fl4.8| ) at rj > 0. According 
to the holographic dictionary, e = l/r(rj) is to be interpreted as a field theory UV cutoff. 
This provides the holographic counterpart to the necessity of introducing a UV cutoff for 
defining the entanglement entropy [39]. Hence 

C„ , c f 1 ' 2 , , 2 



S(l, t) = -L(l, t) = — r(x) 2 dx . (4.10) 

O <-> r * J-q 

The value of the central charge can be read from Q, c = 3/2G^ } . 

We shall present now the results from the numerical integration of the geodesic equa- 
tions. For convenience we will plot the finite, cutoff independent contribution to the 
geodesic length defined by 

2 f 1 / 2 



2 r /z 

L(l,t) = — r(x) 2 dx - 21og r(r/) . (4.11) 
r * Jr] 
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Figure 3: L(l,t) as a function of I for t — 0, 1, ..,8 from bottom to top, for a = 1/3 (left) and 
a = 2 (right). The dashed line corresponds to thermal equilibrium and the dot-dashed line to the 
vacuum result. 



In Fig.|| we show this quantity as a function of the interval size at fixed time, for several 
times that cover from the beginning of the perturbation until well after it has ceased. The 
logarithmic dependence proper of the CFT vacuum (2.1) is reproduced for small intervals, 
since the associated geodesies stay close to the boundary and perceive an AdS geometry. 
A regime of linear growth is established at intermediate values of I some time after the 
action on the system ends. The slope in this regime is proportional to the final temperature 
characterizing the black hole formation process. However for perturbations giving rise to 
a slowly varying geometry, as in FigJ^b, the geodesic length is able to reflect the building 
up of m{v). For large intervals there is an asymptotic return to a logarithmic dependence 
on I. The effect of the perturbation is encoded in a time dependent shift 

S(l,t) -►*(<) + £ log-. (4.12) 
3 e 

The transition from the linear to the logarithmic regime happens at I ~ 2t, being sharper 
the shorter the time extent of the perturbation is. This is clearly seen comparing Fig.^a 
and |b. 

We will denote by v* the value of v attained at the geodesic midpoint. Its dependence 
on the size of the interval is plotted in Fig.[|a. As advanced in the previous section, v* 
decreases as I increases. For the case we are considering of a non-compact x direction, 
v* — > — oo when the interval tends to cover the complete system. Hence, no matter how 
late a field theory time we consider, the geodesies (4.4) are able to probe the whole collapse 
process. 

The linear regime of S(l, t) as a function of I can be easily understood. Although the 
mass function (4.2) has non-compact support, it reaches 98% of its final value at v ~ 2a. 
From then on it is reasonable to consider the black hole as already formed and described by 
a BTZ solution to a very good approximation. BTZ geodesies with the boundary condition 
(4.4) satisfy the very simple relation I = 2(t — i>*). 2 Therefore the condition for geodesies 



This can be deduced from 



v = t - -Z- ln r + ^, i = J=ln r - + V g. (4.13) 
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Figure 4: For a — 1/3 and t — 0,2, .., 10: (left) value of v at the geodesic midpoint as a function 
of the interval size, and on the inset its derivative div*; (right) b in (4.9) as a function of I (solid 
lines) and its asymptotic value 6oo (dashed lines), whose evolution is plotted on the inset. 



to lie entirely in the BTZ part of the geometry is 

l<2t-Aa. (4.14) 

Those that reach the black hole horizon extend tangent to it, and induce a linear depen- 
dence of S(l, t) with the interval size. The BTZ geodesies reproduce the result for the 
entanglement entropy at thermal equilibrium ( |2.2p [17,18]. From that expression we ob- 
serve that the linear behavior requires intervals of size I > /?/2, where (3 = 2tt / \fm is the 
inverse temperature. For them 



Compatibility with the bound (4.14) implies that a clear extensive regime for the cutoff 



independent piece of the entanglement entropy is not established before 

t ~f3/4 + 2a. (4.16) 

The influence of the parameter a on to can be checked in Fig.||. 

Geodesies associated to intervals larger than ( |4. 14 ) enter the black hole formation 



region and present qualitative new features. In this region the apparent horizon builds up 
and clearly differs from the event horizon. Let us describe first the behavior of geodesies 
with endpoints at i>io- When their central segment starts feeling the varying geometry, 
it develops a dip that however does not reach the apparent horizon. This region is followed 
by a segment that crosses behind the apparent horizon and extends at almost constant r 
into the BTZ region, see FigJ||. Both segments lie behind the event horizon. This behavior 
becomes more pronounced as the size of the interval increases. Once v* reaches the pure 
AdS part of the geometry, the dip around the midpoint broadens. The transition to large 
values of the radial coordinate is now through a region well behind the apparent horizon, 



The fist equation is derived from (3.9), and substituting it into the integral of motion (4.6) with m(v)—m 



the second is obtained after integration. 
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Figure 5: Profile of geodesies for different I at i = 6 and a = 1/3. Projection into the (a) r — x 
plane, (b) r — v plane and (c) v — x plane, (d) Full profile of the geodesies. The dashed line in (b) 
and the shaded surface in (d) signal the location of the apparent horizon. In (a) the dashed lines 
represent the position of the apparent horizon at the value of v defined by the trajectory of each 
geodesic: \J m(v(x)). 



where a local minimum in r develops. It is interesting to point out that geodesies that enter 
the black hole formation region require the integration constant b in fl4.9D to be different 
from zero, as shown in Fig.||b. 3 Hence, although a segment of the geodesic lives in the 
BTZ part of the geometry, it differs from that of an equal time BTZ geodesic. This is the 
technical reason that allows these geodesies to cross behind the horizon already in the BTZ 
region. 

For asymptotically large intervals the AdS segment tends to cover the whole extent 
of the geodesic in the x coordinate, explaining the return to a logarithmic growth of the 
entanglement entropy. The piece behind the apparent horizon lies then mostly in the (r, v)- 
plane as can be seen in Figj^b. Indeed, the projection of this piece in the x-coordinates 



The dependence of the bound (4.14) on the parameter a can be neatly seen in this figure, since b starts 
differing from zero before 2t. 
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Figure 6: (left) Radial projection of the infinite interval geodesies as a function of the rescaled 
variable x for a= 1/3 and t = 3, 3.5, 4, 4.5, 5 from left to right. The dashed line represents the posi- 
tion of the apparent horizon as defined in Fig.||a, which practically coincide for all these geodesies. 
The dot-dashed line corresponds to the infinite AdS3 geodesic, (right) Full profile of several t = 6 
long geodesies. 



localizes at x = 0, 1 with extent proportional to r*. This is precisely what is needed for this 
segment to contribute a finite amount to the length functional (4.8) in the limit r* — >0, and 
account for the function s(t) in ( |4.12j ). The limiting large I behavior can be best studied 
by rescaling x — > x = x/r* and setting r* = in equations (|4.6|), fl4.7| ) 



r 4 = 2r'v' — (r 2 — m(v))v 



12 



2vY 



2 12 

r v 



rv 



0. 



(4.17) 



These equations should be integrated with the boundary conditions r(0) = oo, v(0) = t. 
The remaining integration constant is fixed by requiring that the solution reproduces a 
fixed time geodesic in AdS3 far from the boundary, see Fig.^a. It is natural to identify the 
length of this AdS3 geodesic with the I— too limit of the logarithmic term in ( 4.12 ). The 
function s(t) can then be represented in terms of the length difference between the Vaidya 
and the AdS3 limiting geodesies 



s(t) = -AL^i) = £ 

O 6 



r{xfdx-\n r -^$- 



(4.18) 



The point £ has been chosen in the region where both geodesies coincide. The AdS3 limiting 
geodesic is described by 

i i 

(4.19) 



1 

r 



where the matching to the Vaidya geodesic far from the boundary determines the endpoint 
xq. The second term of the bracket in ( |4 . 1 8| ) gives the regularized length of the AdS3 
geodesic segment between £ and the boundary. 

To complete our analysis, we consider the early time behavior of the geodesies. Con- 
trary to what happens for t>to, at early times the geodesies enter the black hole formation 
region without necessarily reaching the apparent horizon. Moreover there is a threshold 



-14- 




A 



Figure 7: (left) Infinite interval geodesies at t = 0.8, 1, 1.4, 1.7 and a = 1/3 from left to right and 
the associated projection of the apparent horizon as in Fig.|]a from right to left (dashed lines). The 
dot-dashed line corresponds to the infinite AdS3 geodesic, (right) Projection in the r — x plane of 
several geodesies for £=1.4. 

time before which geodesies stay always outside it. This happens after the onset of the per- 
turbation. Taking as an example (3 = 2tt and a = 1/3, from Fig we observe that geodesies 
of any length stop crossing the apparent horizon before £~1.4, well after the perturbation 
has ended. At that time S(l,t) already differs appreciably from the vacuum result, with 
s(t) being of order c/6 as can be seen in Fig.|||a. 

Although the apparent horizon is not a necessary ingredient in the holographic deriva- 
tion of the entanglement entropy, its role in it could indicate how much information it 
encodes about the entanglement properties of the field theory. An important related ques- 
tion is wether the apparent horizon can be of relevance to define a notion of field theory 
entropy in far from equilibrium situations [19,40]. Contrary to the case of the event horizon 
in static backgrounds, we have seen that the apparent horizon does not act as a wall for 
the geodesies involved in the calculation of the entanglement entropy. However for t > to it 
strongly influences the geodesic shape and length, with s(t) measuring the geodesic segment 
behind it. Not even this applies in the early, far from equilibrium stages of the evolution, 
when the apparent horizon does not seem to play a role in what regards the holographic 
derivation of the entanglement entropy. 

5. Propagation of entanglement 

The results of the previous section show a strong analogy between our holographic model 
and the critical quantum quenches studied in [1]. In both cases the cutoff independent 
piece of the entanglement entropy behaves after a certain time as an extensive quantity 
up to intervals of size I ~ 2t, see equation fl2.3| ) and Fig.||. Differences arise however for 
larger intervals. For them S(l, t) turns out to reproduce the entanglement properties of the 
initial sea of the excitations generated by the perturbation that brings the system out of 
equilibrium. The excited state created by a quantum quench has a mass gap. Therefore 
only excitations produced at close points can be entangled, and in that case S(l, t) saturates 
to a constant value for I > 2t. The Vaidya geometry models a perturbation on the CFT 
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vacuum which acts during a time <5i~4a. Remarkably the quantum entanglement among 
the resulting excitations appears to be characterized by the same logarithmic scaling as in 
the CFT vacuum. Indeed at t = 0, when the perturbation is half way of being completed, 
S(l,0) ~ §logZ up to a subleading term of order 0(a 2 ). This implies the presence of 
long range correlations among the excitations sourced by the perturbation, contrary to a 
quantum quench. However as in that case we have seen that for all times the entanglement 
entropy of large intervals reproduces the initial entanglement pattern, which now implies 
the return to a logarithmic growth. 

The preservation of the logarithmic scaling of the entanglement entropy during the 
perturbation, and afterwards for large intervals, suggests that the complete evolution in 
the holographic equilibration process is unitary, as it happens in a quantum quench. Since 
the field theory starts in a pure state, a unitary evolution will preserve this character. We 
can use the property that in a pure state the entanglement entropy of a region is equal of 
that of its complementary, as stated in (|1.2| ), to study this question. According to [41], 
the geodesic whose length should give the entanglement entropy of a certain region does 
not need to be connected. It is only required that this geodesic, together with the region 
under consideration, defines the boundary of a certain surface in the bulk geometry. For 
the infinite, complementary region to a finite interval we have clearly two choices [42]. A 
single geodesic with the same boundary conditions (|4.4j ) but extending towards x—> ±00. 
The second possibility is composed of two disjoint pieces: the finite size geodesic whose 
length gives S(l,t), plus a solution of the equations of motion that covers x 6 (—00,00) 
without reaching the AdS boundary. Such a solution is provided by r = at v = — 00, where 
the mass function (|4.2|) strictly vanishes and hence there is no curvature singularity at the 



origin of the radial coordinate. Substituting in (4.5) we obtain that it has zero length, 
as it would be the case in empty AdS. Hence the disconnected choice of geodesies gives 
the minimal length and is the preferred one. It implies that the entanglement entropy of 
a finite interval equals that of its complementary region for any time, proving that the 
evolution is unitary. 4 

The evolution of the entanglement entropy for several fixed intervals is plotted in 
Fig.||. We observe that S(l,t) for a given I attains its final value within a finite time. The 
holographic explanation for this fact is that after a finite time the associated geodesies lie in 
the BTZ part of the geometry and thus produce a time- independent result. An estimation 



for the time at which this happens is obtained from ( 4.14| ) 

t~ l - + 2a. (5.1) 

Except for the dependence on the parameter a, whose effect is clearly noticed in Fig.|8|, 
the same result holds for critical quantum quenches. According to the picture in Fig.|l], 
both the extensive behavior of the entanglement entropy for I < 2t and its saturation as a 
function of time for t>l/2 lead to the important conclusion that entanglement propagates 



4 In the previous section we have seen that geodesies associated to infinite size intervals anchoring in 
the AdS boundary at arbitrary finite time reach v = — 00, as shown in Fig.^i. The same happens in static 
backgrounds. This offers additional support for the previous construction. 
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L(l,t) 



L(l,t) 




Figure 8: L{l,t) as a function of £ for 1 — 2, .., 18, for a = l/3 (left) and a = 2 (right) from bottom 
to top. 



in 2-dimensional CFT's with velocity v 2 = 1 [1]. The Vaidya model allows to check this 
property in an equilibration process with initial long range entanglement. 

It is also interesting to analyze the time evolution of the function s(t). Equation (|4.1Sj ) 
allows its precise numerical evaluation for any time. In Fig.||a we have plotted s(t) for 
several perturbations with different time extent. Its evolution is dominated by a linear 
growth with twice the slope characterizing the extensive regime ( [4.15| ) of the entanglement 
entropy. The linear behavior and the approach to it are well described by 



■<t) = — (2t + ae— 



so 



(5.2) 



where so, ot and 7 are constants. 





afT 



Figure 9: (left) ALoo(t) in fl4.18|) for a = 1/3, 1,2,4, and on the inset its time derivative. As in 
previous figures, f3 = 2ir. (right) The dimensionless quotients a/3 -1 and 7/? _1 as a function of a/3 -1 . 



The dependence of a and 7 on the time extent of the perturbation has been plotted 
in Fig.|9|b. When a<0.2 /3, both a and 7 are independent of a, implying the same property 
for the time needed to achieve the linear growth of s(t). Requiring the term in brackets on 
the rhs of ( |5.2|) to behave linearly within 1% error leads to a threshold time t = r y+Aa. For 
short perturbations we obtain 

f~ 0.85/3, (5.3) 
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The fact that ( |5.3[ ) is independent of the time extent of the perturbation suggests that it 
has a meaning as a characteristic time scale for the evolution of the holographic plasma. 
Notice that t occurs latter than to, the time at which the cutoff independent piece of the 
entanglement entropy starts exhibiting an extensive regime. In the next section we will 
discuss a possible interpretation for (|5.3|). For long perturbations a > 0.5/3, we have instead 
a = a/2 while 7 tends to a constant. Hence in this case the threshold time for the linear 
growth of s(t) is governed by the time extent of the perturbation, and i~io- 

After i not only the asymptotic quantity s(t), but the entanglement entropy itself 
simplifies. Using ( 4.12 ), the following relation is always fulfilled for asymptotically large 
intervals 

(2d l + d t )(S(l,t)-s(t)) = 0. (5.4) 

For t > i this equation extends its region of validity to intervals as small as (3/2, the 
threshold for the onset of the extensive regime, as can be seen in Fig. 10 a,. Integrating we 
obtain 

S(l,t) = f{l-2t) + s{t), t>i,l>0/2, (5.5) 

which applies to the interesting region I ~ 2t where the entanglement entropy changes 
from the extensive to the logarithmic behavior. It implies that after i a regime of steady 
evolution has been achieved, in which the only change of S(l,t) around / ~ 2t is a linear 
increase with time described by the function s(t). Reproducing the extensive regime of the 
entanglement entropy (fi.15 ) determines 

C f3 7TCX 

/( x ) = -log — + — + s , x<-4a, (5.6) 
At large positive values of x the function f(x) inherits from ( |4.12| ) a logarithmic behavior 

f(x) ~ ^ log X - X ° . (5.7) 
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As an example, for f3 = 2ir and a = 1/3 we have obtained xo ~ 7 which provides a good 
description already for x>15. 



Equations (5.2) and ( p75| ) also apply to the evolution after a quantum quench (P^q), 
with the only difference that /(x) = for x > while ( |5.6[ ) holds for x < 0. The analysis 
of quantum quenches performed in [1] requires t, I 3> tq ~ /3 e // an d hence is valid for 
times well after the analogue of t and intervals over the threshold for extensive behavior. 
This is consistent with the regime of validity of (5.2). In this way we found not only a 



qualitative but also quantitative agreement between the evolution of entanglement entropy 
in critical quantum quenches and in our holographic model. The fact that they correspond 
to quite different type of perturbations strongly supports that a very generic behavior for 
the entanglement entropy of 2-dimensional CFT's holds even in time dependent situations. 

6. Comments on thermalization 

A very interesting aspect of the AdS/CFT correspondence is that it allows to model pertur- 
bations which bring the boundary theory far from equilibrium. As for quantum quenches, 
this process and the subsequent evolution are unitary. Hence if the system starts in a pure 
state it will remain so for all times, implying that information cannot be lost at the micro- 
scopical level. However information will become ever more difficult to access as a result of 
the propagation and interaction among the excitations sourced by the initial perturbation. 
This process will drive the system towards a state that at the macroscopical level cannot be 
distinguished from that of thermal equilibrium. In this section we want to further discuss 
what can be learnt about the equilibration process and the holographic dictionary implied, 
from the study of the entanglement entropy. 

When a system is in a pure state, the entanglement entropy gives a measure of the 
quantum entanglement between the degrees of freedom inside and outside a chosen sub- 
system [22-24,26]. In the case of a dynamical situation, it describes also its evolution. 
Propagation will make excitations initially localized in a small region spread over ever 
larger scales. Hence the entanglement entropy of a subsystem will strongly depend on the 
relation between its size and the typical separation of entangled degrees of freedom at a 
given time. This effect is very neatly seen in 2-dimensional CFT's, where the study of 
quantum quenches indicates that this separation evolves as a light front [1]. This was also 
confirmed by our holographic model, whose main difference with a quantum quench is the 
presence of long range entanglement among the initial sea of excitations. 

Interaction is responsible for the redistribution of energy between the different momen- 
tum modes, and therefore it is crucial for the system to relax towards a state of thermal 
equilibrium. Let us consider the entanglement entropy of a 2-dimensional CFT at thermal 
equilibrium (|2.2| ) and perform its Fourier transform. We obtain 



2 jTc=.(pow|— | 



2 eP\p\ - 1 



+ |log^W. (6.1) 



Remarkably the term on the rhs inside the square bracket reproduces the equal time Green's 
function of a free massless scalar at thermal equilibrium. Indeed, the second term in the 
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parenthesis gives the Bose-Einstein statistics at inverse temperature f3. We consider this 
an indication that the entanglement entropy is sensitive to the mode distribution of energy. 
The dispersion relation e(p) = \p\ in ( |6.lD could be directly related to the velocity v 2 = 1 
with which the entanglement appears to propagate in a 2-dimensional CFT. 

Up to very suppressed effects, the action that brings the system out of equilibrium in 
the Vaidya model takes place for t E [—2a, 2a]. If entanglement propagates towards large 
scales with v 2 = l, all entangled excitations at time t will be separated at least a distance 
2t—Aa. This bound has also a holographic meaning: we have shown in ( 4.1 4| ) that geodesies 



associated to smaller intervals only perceive the black hole already formed. For t>to, with 
to in ( |4.16| ), some of these geodesies are able to reach the black hole horizon and induce 
an extensive behavior for the cutoff independent piece of the entanglement entropy. This 
suggests that after to regions of size I < 2t — 4a have locally relaxed to a stationary state 
that can be identified as thermal equilibrium. Of course, in order to better sustain this 
statement it will be necessary the evaluation of other observables. 

Larger intervals will contain entangled excitations. The holographic dictionary for the 
Vaidya model relates such intervals with geodesies that probe the black hole formation 
region, and for I sufficiently large the empty AdS geometry previous to the collapse. In- 
terestingly these geodesies present features absent in their BTZ counterpart. They cross 
behind the event horizon and moreover also behind the apparent horizon. The entangle- 
ment pattern among the sea of excitations and its causal propagation, determines that 
observables involving a length scale larger than 2t — 4a at time t can have an expectation 
value far from thermal. Hence equilibration cannot be achieved at the global level in a 
system of infinite size [1,43]. We relate this property with the necessity to reach behind 
the apparent horizon in the holographic description of the field theory evolution. This 
behavior seems very generic for evolutions which admit a dual representation in terms of 
a process of gravitational collapse. Another observable able to feel that equilibration does 
not reach to regions of size I > 2t are two-point functions, whose evolution after a quantum 
quench was studied in [27,28]. It would be very interesting to explore the behavior of 
two-point functions in dynamical holographic models. 5 

Although we do not know the initial mode distribution created by the perturbation, 
the fact that the Vaidya metric only differs from AdS deep in the interior together with the 
absence of other characteristic scales in the problem, suggests that the typical momentum 
of the excitations sourced by the perturbation is set by the energy density e(t) ([4.3j). This 
would imply that the energy is concentrated on modes of momentum up to T = 1//3, 
with larger momenta being suppressed already at the moment of production. In order 
to describe entangled excitations as separated a certain distance, it is necessary for the 
separation to be bigger than their wavelength. Therefore we can expect that this does not 
clearly happen before to, especially for short perturbations. On the gravity side, geodesies 
for t < to enter the black hole formation region without necessarily crossing the apparent 



5 Based on semiclassical arguments, it has been proposed that the spacelike two-point function of oper- 
ators with very high conformal dimension A is given holographically by e~ AL , with L the length of the 
geodesic that joints the two insertion points [7,44]. For these reduced set of two-point functions, our analysis 
of geodesies implies consistency with the results of [27, 28] . 
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horizon. As t decreases the apparent horizon influences ever less the shape of the geodesies. 
Moreover there is a time before which the geodesies never reach the apparent horizon, 
while the entanglement entropy already deviates appreciably from the vacuum result. For 
perturbations with a short time extent this happens well after the perturbation has ended. 
Hence the apparent horizon seems to lose relevance in the early stages of the evolution, as 
the picture of separated entangled excitations becomes increasingly inaccurate. 



Besides to an important time in the evolution of our holographic plasma is i, after 
which the entanglement entropy enters a regime of steady evolution as described in (|5.5| ) . 
We have argued based on (|6.1|) that the entanglement entropy is sensitive to the mode 



distribution of energy. Assuming this, a possible explanation for the regime (5.5) could be 
that the mode occupation numbers have reached their final distribution. The only dynamics 
remaining on the system would be associated to the propagation of entanglement over ever 
larger scales, a process whose evolution is constrained by causality. For perturbations with 
a short time extent, i is independent of a and of the same order but larger than to, the 
time at which local thermalization can be first identified. Both facts suggest that ([^) is 
a characteristic time for the evolution of the holographic plasma, not directly related to 
the propagation of entanglement. This adds support to the interpretation of i in (5J3) as 
a thermalization time for occupation numbers and would imply a very fast thermalization 
at strong coupling, on the line of related results for holographic plasmas [10-12,45]. At 
any rate, a firmer statement about the evolution of occupation numbers would require 
analyzing better adapted observables. 



Finally we would like to show that a simple model similar to that proposed in [1] to 
reproduce the evolution of entanglement entropy after a quantum quench for Z,i^>0 (|0|), 
can be generalized to account for the Vaidya evolution after t. This model assumes that 
left and right moving entangled excitations propagate classically and without scattering, 
with a velocity v 2 = 1 in the case of a 2-dimensional CFT. Although this hypothesis 
could find support from a quasiparticle picture at weak coupling, such justification would 
not apply to strong coupling. However the evolution ( [2 .3D and therefore the success of 
the model is independent of the strength of the coupling. This suggests to drop any 
quasiparticle interpretation and consider the model as an effective way to describe the 
steady propagation of entanglement. Interestingly, the propagation of entanglement with 
velocity v 2 = 1 is supported not only by the study of dynamical setups, but by the result 
of Fourier transforming the entanglement entropy at thermal equilibrium ( |6.1| ) . 

We generalize the model in [1] by allowing for the presence of entanglement among 
excitations originating from points separated a distance y. For the sake of simplicity 
we will focus on a very brief perturbation, a ~ 0. We introduce a function p(y) that 
describes the density of entangled excitations separated a distance y at the moment of 
their production. Given the translational invariance of the perturbations in our model, 
this function is independent of their location. In order to implement the steady regime 
( |5.5| ), we assume that at time t the entanglement among degrees of freedom at a distance 
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y is described by 



U y < 2t . 



Correlated excitations at x and x + y will contribute to the entanglement entropy of an 
interval when only one of these points lies inside the chosen interval 

pod />oo 

S(l,t)= dx dy(h(x)-h(x + y)) 2 p(y,t), (6.3) 



with h(x) = 1 if x G [0,/] and zero otherwise. Substituting (|6.2|) , we obtain an extensive 
result for intervals of size I < 2t 



oc 



S(l,t) = 2l p{y)dy. (6.4) 



o 



Coincidence with the cutoff independent term of ( |4.15 ) determines pdy = 7rc/6/3. For 
l>2t we then have 

s ( l > *) = -&r + 2 ^ - 2t ) / + / 2 yp(y) d y > ( 6 - 5 ) 

in agreement with ([T^). Moreover from the above relation we derive 

p(y) = -ldfS(l,t)\ l=2t+y . (6.6) 

The return of S(l,t) to a logarithmic growth for sufficiently large intervals in the Vaidya 
model implies indeed the quantum entanglement between excitations produced with arbi- 



trary large separations. Substituting (|5.5[) and (5/7) we obtain 

P(y) = -\d 2 y f(y) -+ I . (6.7) 

A preliminary analysis indicates that xq has only a mild dependence on the parameter a 
and thus is approximately proportional to f3. For the evolution after a quantum quench 
( |2.3j ), relation fl6.6|) leads to p(y) = ^S(y). Therefore we recover the initial localization of 
entanglement on very small scales and the picture in Fig.|l]. 

In this paper we have explored the application of the AdS/CFT correspondence to the 
study of far from equilibrium dynamics in strongly coupled field theories, a very topical 
subject due to the heavy ion experiments at RHIC and LHC. We have focussed on the 
evaluation of the entanglement entropy, an observable with a simple holographic represen- 
tation [17,18] generalizable to time-dependent situations [19]. A natural extension of our 
work would be to study the evolution of entanglement entropy for field theories in dimen- 
sion higher than two, where few results are available. It would be also very interesting 
to evaluate the entanglement entropy in an equilibration process dual to the gravitational 
collapse of a shell with a more natural composition than null dust, as for example it was 
considered in [13]. We plan to address some of the issues and open questions discussed in 
this section in future work. 
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